Exact Thermodynamics of the Double sinh-Gordon Theory in 1+1-Dimensions 
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We study the classical thermodynamics of a 1 + 1-dimensional double- well sinh-Gordon the- 
ory. Remarkably, the Schrodinger-like equation resulting from the transfer integral method 
is quasi-exactly solvable at several temperatures. This allows exact calculation of the parti- 
tion function and some correlation functions above and below the short-range order ( "kink" ) 
transition, in striking agreement with high resolution Langevin simulations. Interesting con- 
nections with the Landau- Ginzburg and double sine-Gordon models are also established. 
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The statistical mechanics of nonlinear coherent struc- 
tures in low dimensions has long attracted theoretical 
attention, both for the intrinsic interest in such funda- 
mental problems as kink nucleation and dynamics, as 
well as in diverse applications, e.g., in conducting poly- 
mer physics 0] and DNA denaturation ||^. Both ana- 
lytic and numerical techniques have been applied to these 
problems: well known among them are the (analytic) 
transfer integral method and the (numerical) Langevin 
method. The transfer integral technique converts the 
problem of finding the classical partition function Zd to 
an eigenvalue problem for a Schrodinger-like equation to 
which familiar approximation methods such as WKB can 
then be applied Q . The advantage of Langevin methods 
is that (unlike Monte Carlo) real time quantities such 
as temporal correlation functions can be computed, and 
kinks/ antikinks tracked both in space and time. 

In the past, applications of these methods have yielded 
comparisons of approximate analytic results with numer- 
ical data to only rather low levels of accuracy, of order 
of tens of per cent Q. In this Letter we report substan- 
tial progress on both fronts. We discuss a nonintegrable 
1 -|- 1-dimensional field theory for which thermodynamic 
quantities can be computed exactly at several tempera- 
tures using techniques from quasi-exactly solvable (QES) 
potentials in quantum mechanics. This theory is in the 
same class as the more familiar Landau-Ginzburg model 
and also admits exactly known kink solutions ||6[. We 
have carried out very high resolution Langevin simula- 
tions and find excellent agreement with the exact results 
at the checkpoint temperatures. The high accuracy of 
the Langevin simulations allows the use of the probabil- 
ity distribution function (PDF) to directly compute ther- 
modynamic quantities m thus providing an alternative 
to conventional methods based on fluctuations. 

In order to calculate Zd exactly at various temper- 
atures one has to solve a Schrodinger equation with a 
temperature dependent mass. While completely solvable 
potentials are rare, in the last few years several double- 
well QES models have been discovered for which the ex- 
act classical partition function can be found at one given 



temperature. The drawback has been that the exact 
eigenstates are only known for a given set of couplings 
and as a result, it has not been possible to obtain the 
exact Zci at more temperatures. In this Letter we show 
that for the double sinh-Gordon (DSHG) QES problem, 
if the ground state energy is known for n different values 
of coupling constants, then Zd can be evaluated for any 
of these theories (with a given set of coupling constants) 
at n different temperatures. This is also true for the 
triple well (f>^ and the double sine-Gordon (DSG) mod- 
els, results for which will be reported elsewhere M. We 
conjecture that this result holds for a large class of QES 
problems. 

The double-well (p'^ model in 1 4- 1-dimensions has 
been extensively studied. However, in this case the 
Schrodinger equation does not possess known exact solu- 
tions. To overcome this problem we turn to the DSHG 
potential: 
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{(/)) = (Ccosh2(/)- n)^ 
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where ^ is a positive parameter. In order to have a 
double-well potential, n > ^, in which case the two min- 
ima are located at cosh20o = Moreover, for the 
system to be QES, n has to be a positive integer. This 
potential is the hyperbolic analog of the double sine- 
Gordon system. Similar potentials arise in the context 
of the quantum theory of molecules {e.g. a homonuclear 
diatomic molecule), wave motion describing the normal 
modes of vibration of a stretched membrane of variable 
density, and as the solution of a Fokker-Planck equation 
Q. The hyperbolic analog of the sine-Gordon equation 
is a single well potential (sinh-Gordon) and thus uninter- 
esting from the soliton statistical mechanics perspective. 

The DSHG potential written in the form (^ has all 
the generic features of a double-well potential such as 
Landau-Ginzburg, but allows for much greater analytic 
progress. Below we find exact solutions for 1) a kink, 
2) phonon dispersion, 3) a kink lattice, and 4) the first 
few eigenvalues and eigenfunctions of the transfer oper- 
ator at certain temperatures, allowing thereby analytic 
calculation of the PDF and correlation functions in the 
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thermodynamic limit. 

We exhibit below the exact kink and kink lattice so- 
lutions for the DSHG theory (details will be given in 
Ref. 1^). A kink is a time independent solution result- 
ing from the minimization of the total energy density 
e{x) — Vdshg{4>) + (.9/2)01 with the boundary condi- 
tions (j) — > ±(/)o as X — > ±cxD. The constant g is often 
introduced in condensed matter treatments as a phe- 
nomenological parameter and controls the kink size. In 
a field theoretic context, however, g = \, and this is the 
value we choose here. (All the solutions given below can 
be written for arbitrary g.) 

The kink/antikink solution, located at xo is, (n > C,) 



{x) = ± tanh tanh 0o tanh 



(2) 



where tanh^o = 0/(^ + 0, f = [2(^-0]-^'"^. 

The kink has topological charge Q = ^dx = 20o- 
Traveling kink solutions are obtained by boosting to ve- 
locity V via X ^ (1 — v'^)~^^'^{x — vt). The statistical 
mechanics of kinks is governed largely by the kink en- 
ergy (or rest mass): 



Eg = 4^n^n2 - (2 tanh 
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The phonon dispersion around the minima ±0o for this 
model is = q"^ + S{n'^ - C^) = + and the 

phonon contribution to the free energy per unit length is 



-t 1)'. 
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with 5 being the lattice constant and /3 = 1/ /cbT. 

In order to understand kink-antikink interactions, it 
is very useful to construct kink lattice solutions (a 
kink/antikink chain). For the DSHG theory, this solu- 
tion is 



,{x) = ± tanh ^ ( tanh 01 sn 
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where d is the periodicity of the kink lattice, K{k) is the 
complete elliptic integral of the first kind with modulus 
fc, sn(x, k) is the Jacobi elliptic function, and with < 
< ^0 = 0) = (n - C)^ 



cosh2(/>i,2 = T 



C 



(5) 



The topological charge (per period) in the lattice problem 
Ql = 2(j>L{K) = 2(j)i is smaller than the single kink case. 
The kink size in the kink lattice, ^l, is also smaller than 
the free kink size 



The energy of the kink lattice per period {i.e. energy 
per kink-antikink pair plus the interaction energy) is 



El = 4a 



f2 

(n + CfK + ^{n^ - C){K - E) - 4nCn 



where E{k) and n(tanh^0i,fc) are complete elliptic in- 
tegrals of second and third kind, respectively. In the 
dilute limit (fc — > 1, d — > 00) the divergences in K{k) and 
n(tanh^ 01, k) exactly cancel out and we recover the sin- 
gle kink result Eg. The interaction energy as a function 
of separation [i.e. k or d) is given by Ein = — 2Es. 

Turning now to the computation of Zd, we note that 
this calculation can be divided into two parts: a triv- 
ial Gaussian integration over the field momentum, and 
a computation of the configurational partition function, 
which via the transfer integral method becomes equiv- 
alent to solving a Schrodinger-like equation The 
Hamiltonian for the DSHG theory is 



H 



dx 



\tt^ + \{d,4>? + Vdshg{4>) 



(6) 



and this leads to the Schrodinger equation for the eigen- 
values and eigenfunctions of the transfer operator. 



-^^v,, + (Ccosh20- 



n)Hk = Ek'^k (7) 



Remarkably, this equation is an example of a QES sys- 
tem. Using results for a related potential from Ref. |^ , at 
2/3^ = 1 the eigenstates of the first n levels can be found 
for n — 1,2,3,4. (We have extended this to the cases 
n = 5,6.) However, what one really wants is to consider 
a given fixed-n theory and obtain eigenstates at different 
temperatures. It is easy to see from Eq. (^, by simple 
rescaling, that solutions of a fixed-n theory at certain val- 
ues of (3 are the same as the solutions of another theory 
(different n and C) at 20^ = 1. Depending on the chosen 
value of n, exact solutions are available at different fixed 
values of j3. Here, we restrict ourselves to one such fam- 
ily (n — 2) which allows the exact computation of the 
first few eigenstates at 8/3^ — m? [m — l,---,6). For 
illustration, two examples of the (unnormalized) ground 
states are given below (see also Fig. 1). The first (high 
temperature, = 1/8) has an eigenfunction with a sin- 
gle peak while the second (lower temperature, 0^ = 1/2) 
has a double peak: 
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*o(0)|/32^i = exp — -Ccosh20 



^o(0)|/32^i = cosh0exp ^— i(cosh20J , (8) 

with corresponding ground state energies, Eq = 1 -f C^/4, 
Eq ^ -2C + 3. The PDF for the field is the square 
of the normalized ground state eigenfunctions. Solutions 
at higher energies and other values of (3 are given in Ref. 
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Once the eigenvalues of the transfer operator are 
known, they can be used to compute the correlation func- 
tions Ci(x) = {(l){0)(j){x)) andCs = {6 (j)^ (0)6(1)'^ {x)), using 

Ci(x) = J2 I (^'^ 1-^1*0)1' exp [-f3\x\{Ek - Eo)] , (9) 

k 

C2{x)=Y,\{'^k\5cf>^\'i'o)\^^^P[-(3\^\iEk-Eo)]. (10) 

k 

It is apparent that at large distances, Ci and C2 are 
dominated by the lowest state with nonvanishing matrix 
elements: the first excited state in the case of Ci and the 
second excited state in the case of C2. Since Eq, Ei, and 
E2 are known at certain temperatures, the large distance 
behavior of these correlation functions can be found ex- 
actly and compared with the results from simulations. 
Static structure factors may also be calculated in much 
the same way. 

At this point, it is important to mention the connection 
between the "quantum" calculations and kink physics. 
In the context of kink statistical mechanics, it is usual 
to introduce a phenomenological description of kinks as 
particles in a grand canonical ensemble. However, this is 
unnecessary, and all such thermodynamical information 
can be extracted directly from the Schrodinger descrip- 
tion of the transfer operator. For example, the kink den- 
sity has been obtained in this way in Rcf. []l0| . Simpler 
quantities like Ci and C2 have obvious natural interpre- 
tations in terms of kinks. The Ci correlation length is 
related to the kink/antikink spacing and increases mono- 
tonically as P increases (Fig. 2). 

The behavior of the C2 correlation length requires a 
little explanation, since C2 is not directly sensitive to do- 
main size. At both high temperatures (no kinks) and 
low temperatures (number of kinks exponentially sup- 
pressed), the correlation length is essentially that set by 
thermal phonons, and is therefore small. However, at 
temperatures close to the kink transition, nonlinear fluc- 
tuations on the kink length scale become important and 
can dominate C2. At these temperatures one expects 
the C2 correlation length to rise to a maximum value, 
of order the kink size, and this is indeed what we ob- 
serve numeric ally (Fig. 3). The Schottky anomaly in the 
specific heat [|o| arises for the very same reason. 

The exact results described above can be compared 
against those obtained from Langevin methods. The ad- 
ditive noise Langevin equation for the DSHG theory is 

92 = dl^4> - ridt4> - 4C(Ccosh20 - 2) sinh20 + F{x, t), 

where F{x,t) is a stochastic (Gaussian, white) external 
force which satisfies the fluctuation-dissipation relation 
linking the noise strength to the viscosity 77, 

{F{x, t)F{x\ t')) = 2'q(3-^5{x - x')5{t - t'). (11) 

This stochastic PDF can be solved by standard methods 
pT| which we implemented on massively parallel com- 
puters. Typical choices for the lattice constant are often 



dictated by memory limitations rather than by accuracy. 
Comparisons with the exact results have led us to con- 
clude that errors up to 30% may be expected if lattice 
discretization is done as coarsely as has been the norm so 
far in numerical calculations. The existence of nontriv- 
ial exact continuum results has proven to be essential in 
carefully estimating error and convergence in field theo- 
retic Langevin simulations |]l2| . Our present simulations 
were typically performed on 5 x 10^ site lattices with a 
lattice constant S = 0.025 and time-step e = .005. 

Fig. 1 shows the striking agreement between the nu- 
merically obtained and the exact continuum PDFs at 
three temperatures: The worst case departure is at the 
level of parts per thousand. The comparisons for the 
inverse correlation lengths are given in Figs. 2 and 3. 
DSHG system parameters are n = 2, — 0.05. For 
Ci, the numerical values are 1/A = 0.1425 (/3^ = 1/2) 
and 1/A = 0.012 (/3^ = 9/8) as compared to the exact 
values in the continuum theory of 0.14142 and .0105, re- 
spectively. The small offset between the continuum and 
lattice calculations is due to the finite value of the lattice 
constant and is consistent with estimates from higher- 
order contributions to the transfer integral ]T^. 
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FIG. 1. The numerically evaluated PDFs at three values of 
/3^: 1/8 (triangles), 1/2 (diamonds), and 9/8 (squares). The 
corresponding continuum exact solutions are the solid lines. 
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FIG. 2. The numerically obtained inverse correlation 
lengths from Ci{x) for fJ'^: 1/2 (diamonds), and 9/8 (squares). 
The large |a;| continuum exact results are the solid lines. 
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FIG. 3. The numerically obtained inverse correlation 
lengths from C2 (2;) for the same three temperatures and with 
the same conventions as in Fig. 1. The solid line is the large 
continuum exact result for 0^ = 9/8. The largest correla- 
tion length is at the intermediate value of jS (see text). 

The high quality of these numerical simulations im- 
plies that the PDF can now be used directly to compute 
thermodynamic quantities at any temperature. Since the 
PDF is just the square of the ground state wave func- 
tion of the Schrodinger equation one can use it to 
compute the ground state energy Eq numerically, from 
which the internal energy {U = dEg/dP), the free energy 
(F = Eo/P), and the entropy {S = PdEoJd(3 - Eq) can 
all be computed in a straightforward way [pj . The specific 
heat involves two /3 derivatives and is difHcult to obtain 
with good accuracy but in this case, the standard energy 
fluctuation method is quite effective. The use of the PDF 
complements traditional techniques utilizing energy fluc- 
tuations in Langevin simulations which are not suited to 
free energy and entropy calculations. 

The QES nature of the DSHG theory allows not only 
the exact computation of Eq at several temperatures, but 
also of dEo/dp, using flrst order perturbation theory: 
dEo/dp\f3=f3o = (^-0, 9^*0/902) where /3o is one of the 
special temperatures e.g., Eq. (|^). Thus the internal 
energy U and the entropy S can also be found exactly at 
these temperatures Once again, these quantities can 
be used to validate numerical work over a broad range of 
temperatures. 

As a final point, we consider the relationship of the 
DSHG theory to the more familiar Landau-Ginzburg 
model. Scrutiny of Eqs. reveals the following 

important connection between the kink (and kink lat- 
tice) solutions of the 0"* model and the double sine- 
Gordon (DSG) and DSHG models. Consider the 0"* po- 
tential Vi{u) — [{n + C)w^ — {n — C)]^- The substitution 
u = tanh0 takes the (static) equations of motion over 
to the DSHG equations. The alternative substitution 
u — tan (p leads to the DSG model. This means that all 
known solutions of the 0** theory can be directly taken 
over to the DSHG and DSG theories (and vice versa). 
As one use of this interesting relationship, the DSG kink 
lattice solution (not known heretofore in the literature) 
can be written down directly for Vdsg = (C cos 2(/)— n)^: 



(f>L—±tan ^ ^tan^i sn ^ — - — -,kjj , (12) 

simply by using the substitution tanh tan in Eq. (Q). 

This connection enables us to write down by inspection 
not just the kink solutions but their total energy as well, 
which is often a very tedious task. Moreover, since we 
know that the DSHG model is an example of a QES sys- 
tem, and considering the very similar way in which the 
DSG and DSHG models are related to (j)'^, it is logical to 
conjecture that the DSG model must also be a QES sys- 
tem. Indeed, this is the case, and we have found several 
exact eigenvalues and eigenf unctions for many temper- 
atures. The exact statistical mechanical results for the 
DSG model, similar to the DSHG results presented here, 
will be reported later . 
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